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This  work  highlights  the  gains  of  an  innovative  model-based  control  approach  applied  to  a  proton 
exchange  membrane  fuel  cell  (PEMFC)  system,  included  in  a  stand-alone  hybrid  generator.  This  approach 
proposes  a  multivariable  setpoint  tracking  of  the  PEMFC  output  power  and  temperature.  The  freshness  of 
this  approach  is  based  on  the  combination  of  a  nonlinear  model-based  predictive  control  strategy  (NMPC) 
and  a  global  linearizing  control  (GLC)  algorithm. 

The  performance  of  the  proposed  control  strategy  is  confirmed  thanks  to  simulations  of  varying  control 
scenarios.  Results  show  good  performance  in  setpoint  tracking,  disturbances  rejection  and  robustness 
against  plant/model  mismatch  in  presence  of  noisy  signals.  Moreover,  for  similar  setpoint  point  tracking 
accuracy,  the  proposed  control  strategy  appears  to  be  four  times  faster  than  a  classical  multivariable 
NMPC  strategy.  According  to  real-time  application  objectives,  this  control  strategy  appears  as  a  promising 
option  to  be  implanted  in  the  overall  control  scheme  of  the  stand-alone  hybrid  power  generator. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nowadays,  in  a  context  of  sustainable  development,  and  to  deal 
with  the  serious  problem  of  global  warming,  renewable  energy 
sources  such  as  solar,  wind  or  biomass  resources  represent  a 
promising  alternative  to  the  traditional  fossil  fuels.  In  a  general 
point  of  view  subtropical  region,  solar  radiation  is  high,  which  is 
favorable  for  various  kinds  of  solar  energy  implementations.  In  par¬ 
ticular  case  of  the  tropical  islands  photovoltaic  appears  as  one  best 
strategy  in  their  energy  mix.  Indeed,  photovoltaic  (PV)  technology 
directly  converts  solar  energy  into  electricity,  without  greenhouse 
gases  and  air  pollutant  emissions.  Nevertheless,  this  technology  has 
two  main  drawbacks:  its  low  efficiency  and  the  nature  of  solar 
power,  intermittent  and  unreliable.  In  practice,  considering  the 
zero  cost  of  solar  energy,  the  weak  conversion  rate  is  usually  not 
considered  as  a  major  disadvantage.  However,  PV  power  conver¬ 
sion  is  limited  to  daylight  period  and  strongly  influenced  by  the 
weather. 

To  overcome  these  limitations,  one  of  the  most  widespread 
methods  is  to  integrate  PV  generators  with  other  power  sources 
such  as  fuel  cells,  batteries  backup  or  diesel  generators  [1,2].  Due 
to  its  high  efficiency,  fast  power  response,  high  power  density  and 


*  Corresponding  author.  Tel.:  +33  476826586;  fax:  +33  476826777. 
E-mail  address:  jonathan.deseure@lepmi.grenoble-inp.fr  (J.  Deseure). 

0378-7753 /$  -  see  front  matter  ©  2012  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 2.01. 096 


low  operating  temperature,  the  proton  exchange  membrane  fuel 
cell  (PEMFC)  appears  as  a  suitable  and  reliable  option  to  be  com¬ 
bined  with  photovoltaic  arrays.  Recent  studies  showed  its  great 
potential  to  solve  the  inherent  problem  of  intermittent  power  gen¬ 
eration  of  PV  arrays  [3-5].  Regarding  photovoltaic-fuel  cell  (PVFC) 
hybrid  systems,  the  main  part  of  the  process  control  remains  the 
control  of  the  fuel  cell. 

A  fuel  cell  is  a  nonlinear  and  strongly  coupled  dynamic  system.  It 
is  a  multi-input  multi-output  (MIMO)  system  based  on  multiphase 
flow,  electrochemical  reactions  and  heat  transfer.  Since  few  years, 
many  control  strategies  have  been  developed  to  control  PEMFC 
and  SOFC  (solid  oxide  fuel  cell)  systems,  based  on  various  con¬ 
trol  algorithms,  manipulated  and  controlled  variables.  Huang  et  al. 
[6]  designed  a  fuzzy  PID  controller  to  regulate  hydrogen  flow  rate, 
resulting  in  an  optimization  of  the  hydrogen  consumption.  Wang 
and  Ko  [7]  proposed  a  multivariable  robust  PID  controller,  based  on 
a  classical  PID  approach  combined  to  the  merits  of  robust  control,  to 
regulate  air  and  hydrogen  flow  rates.  Results  obtained  showed  good 
performance  in  terms  of  stability  and  efficiency  using  a  simple  con¬ 
trol  structure.  Methekar  et  al.  [8]  considered  a  MIMO  system  and 
proposed  two  PID  control  strategies.  Using  a  steady-state  relative 
gain  array  (RGA)  analysis,  these  authors  showed  that  hydrogen  and 
coolant  inlet  flow  rates  are  suitable  manipulated  variables  to  con¬ 
trol  the  output  power  and  temperature  of  a  PEMFC,  respectively. 
Zhang  et  al.  [9]  proposed  a  nonlinear  predictive  control  approach 
to  control  output  power,  fuel  utilization  and  temperature  of  a  SOFC 
using  the  current  density,  the  fuel  and  air  flow  rates  are  used  as 
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Nomenclature 

Afc 

cell  electrode  area,  cm2 

Q 1 

double  layer  capacitance  of  single  cell,  F 

Ch2 

hydrogen  concentration  at  the  anode/membrane 
interface,  mol  cm-3 

C02 

oxygen  concentration  at  the  cathode/membrane 
interface,  mol  cm-3 

Cpw 

water  heat  capacity,  J  kg-1  K-1 

ct 

thermal  capacitance,  JK_1 

Ecool 

heat  flux  removed  by  the  heat  exchanger,  W 

Eelec 

power  consumption  by  the  electrical  load,  W 

Eloos 

heat  flux  dissipated  to  the  surroundings,  W 

E Nernst 

thermodynamic  potential,  V 

Etot 

total  power  delivered  by  the  stack,  W 

F 

Faraday  constant,  96485  C  mol-1 

Eamb 

convective  heat  transfer  coefficient,  W  m-2  K-1 

Econd 

heat  exchanger  conduction  index,  WI<-1 

hconv 

heat  exchanger  convection  index,  WK-1  A-1 

I 

cell/stack  current,  A 

J 

cost  function  of  regulator 

k0 . kr 

_  1  controller  tuning  parameter 

Lrfh 

rth  Lie  derivatives  of  h(x(t))  along/ 

Ugh 

rth  Lie  derivatives  of  h(x(t))  along  g. 

lm 

membrane  thickness,  cm 

iflcw 

cooling  water  mass  flow  rate,  kg  s-1 

2, in 

hydrogen  inlet  flow  rate,  mol  s-1 

2, out 

hydrogen  outlet  flow  rate,  mol  s-1 

mU2, reacted  rate  of  hydrogen  consumption,  mol  s-1 

™02,in 

oxygen  inlet  flow  rate,  mol  s-1 

Wo  2,  out 

oxygen  outlet  flow  rate,  mol  s-1 

m0  reacted  rate  of  oxygen  consumption,  mol  s-1 

Ncell 

number  of  cells 

Nu 

control  horizon,  s 

Ny 

predictive  horizon,  s 

Ni 

optimization  parameter 

Ph2 

hydrogen  partial  pressure,  atm 

PL] 

hydrogen  partial  pressure  at  the  inlet,  atm 

Po2 

oxygen  partial  pressure,  atm 

pout 

oxygen  partial  pressure  at  the  outlet,  atm 

r 

relative  degree  of  the  system 

R 

universal  gas  constant,  8.314  J  mol-1 1<-1 

R-act 

activation  resistance,  £2 

S 

external  stack  surface,  m2 

t 

time,  s 

T 

stack  temperature,  K 

Eamb 

ambient  temperature,  I< 

Ec,in 

cooling  water  inlet  temperature,  K 

Tc,out 

cooling  water  outlet  temperature,  K 

U 

input  vector 

u(k) 

control  variable  at  time  k 

Ua 

overall  heat  transfer  coefficient 

Van 

anode  volume,  m3 

Vact 

absolute  value  of  activation  overvoltage,  V 

Vca 

cathode  volume,  m3 

Vcell 

single  cell  voltage,  V 

V stack 

stack  voltage,  V 

Vw 

cooling  water  volume,  m3 

w 

setpoint 

X 

state  vector 

y 

output  vector 

ym 

model  output 

ysp 

setpoint  to  track 

ySr 

target  value 

Greek  letters 

a  filtering  parameter 

a0  controller  tuning  parameter 

P  weighting  parameter 

AH  hydrogen  consumption  enthalpy,  kj  mol-1 

y  weighting  parameter 

r\act  activation  overvoltage,  V 

rj ohm  ohmic  overvoltage,  V 

pw  water  density,  kg  nrr3 

.  .,£4  activation  overvoltage  parametric  coefficients,  V  or 

VK"1 


manipulated  variables.  In  this  approach,  the  exit  gases  composition 
and  temperature  are  estimated  using  moving  horizon  estimation 
(MHE)  methods.  Wu  et  al.  [  1 0]  designed  a  multi-loop  nonlinear  pre¬ 
dictive  control  scheme  to  regulate  stack  temperature  and  oxygen 
excess  ratio  of  a  hybrid  energy  system  (wind  turbine  and  PEMFC) 
by  manipulating  coolant  and  air  mass  flow  rates.  Even  if  the  sim¬ 
ulation  results  are  encouraging,  the  authors  underline  that  some 
devices  and  design  are  not  taken  into  account  (e.g.  hydrogen  stor¬ 
age,  power  management  and  inverters).  Li  et  al.  [11]  proposed  a 
nonlinear  robust  control  of  a  PEMFC  based  on  a  state  feedback  lin¬ 
earizing  approach.  This  well-known  nonlinear  approach  allows  to 
obtain  a  nonlinear  control  law  directly  from  the  dynamic  nonlinear 
model  of  the  process.  To  improve  the  robustness  of  the  system,  the 
authors  proposed  to  add  a  robust  control  strategy  to  the  state 
feedback  control  law. 

This  study  focuses  on  the  control  of  a  5  kW  PEMFC,  included  in  a 
stand-alone  hybrid  renewable  power  generation  system  designed 
for  a  typical  household  in  La  Reunion  Island  located  at  following 
geographic  coordinate  (21°  6'  S/55°  36'  E).  In  this  purpose,  a  MIMO 
dynamic  nonlinear  model  of  the  PEMFC,  dedicated  to  nonlinear 
model-based  control  approaches,  is  proposed.  Directly  based  on  the 
nonlinear  dynamic  model  of  the  system,  a  nonlinear  model-based 
predictive  control  strategy  (NMPC)  and  a  global  linearizing  con¬ 
trol  (GLC)  approach  are  proposed  to  control  the  output  power  and 
the  temperature  of  the  PEMFC,  respectively.  This  original  control 
scheme,  designed  for  on-line  applications,  appears  to  be  suitable  to 
cope  with  large  setpoint  changes,  severe  unmeasured  disturbances 
and  plant/model  mismatch. 

The  subsequent  developments  of  the  manuscript  are  organized 
as  follows:  Section  2  is  dedicated  to  the  description  of  the  sys¬ 
tem  and  the  model  design.  The  third  section  deals  with  the  control 
design,  and  details  the  NMPC  of  the  output  power  and  the  GLC  of  the 
temperature  of  the  stack.  In  Section  4,  the  performance  of  the  pro¬ 
posed  control  strategy,  in  terms  of  setpoint  tracking,  disturbances 
rejection  and  robustness  against  plant/model  mismatch,  is  demon¬ 
strated  via  simulation.  Finally,  conclusions  and  prospects  are  drawn 
in  Section  5. 

2.  System  design  and  modeling 

2A.  System  description 

This  work  deals  with  the  modeling  and  control  of  a  Ballard  5  kW 
PEMFC,  represented  in  Fig.  1. 

This  PEMFC  is  including  in  a  stand-alone  hybrid  system.  This 
hybrid  generator  comprises  a  PV  array,  an  electrolyzer,  a  hydro¬ 
gen  tank  and  a  PEMFC.  The  PV  array  is  designed  to  satisfy  the 
power  demand  of  a  typical  household  in  La  Reunion.  In  the  absence 
of  sufficient  solar  radiation,  the  PEMFC  is  designed  to  compen¬ 
sate  for  the  PV  generator,  assuming  the  hydrogen  storage  problem 
solved.  During  the  daylight  period,  the  solar  power  is  converted  into 
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Fig.  1.  Scheme  of  the  PEMFC  system  (with  heat  exchanger  system). 


2.2.  PEMFC  model 

This  section  deals  with  the  dynamic  nonlinear  model  of  a 
PEMFC  dedicated  to  nonlinear  model-based  control  approaches. 
This  is  based  on  a  set  of  generalized  empirical  equations  devel¬ 
oped  for  a  basic  model  that  could  describe  several  types  of  fuel  cells 
[10,13-15].  In  order  to  simplify  the  analysis,  several  assumptions 
are  made  and  are  listed  below. 

-  All  gases  are  modeled  as  ideal  gases. 

-  The  total  pressure  inside  the  stack  is  uniform. 

-  Gas  streams  are  saturated  with  vapor. 

-  Pure  hydrogen  is  supplied  to  the  anode. 

-  The  air,  that  fed  the  cathode,  is  a  mixture  of  nitrogen  and  oxygen 
by  a  ratio  of  21:79. 

-  The  electrolyte  is  fully  saturated  with  water. 

-  The  stack  is  well  designed  in  such  a  way  cells  perform  similarly 
and  can  be  lumped  as  a  stack. 

-  The  temperature  is  homogeneous  throughout  the  stack. 

-  Gases  enthalpy  variations  are  negligible. 


electric  power.  A  fraction  of  PV  power  is  used  to  reach  the  power 
demand  of  user  load  whereas  excess  energy  is  stored  in  hydrogen 
via  electrolysis.  The  fuel  cell  generator  can  be  used  both  at  night 
to  cover  the  energy  demand,  and  during  day  as  backup  genera¬ 
tor  when  solar  power  is  deficient  (e.g.  during  a  period  of  low  solar 
radiation).  Fig.  2  illustrates  a  typical  daily  power  consumption  of  a 
household  located  in  La  Reunion  Island,  and  the  power  delivered 
by  the  PV  array.  Energy  conversion  via  electrolysis  and  fuel  cell  has 
an  efficiency  of  30-35%  [12].  Considering  an  estimated  household 
electricity  consumption  of  112  Whm-2  day-1  and  using  polycrys¬ 
talline  PV  cells  with  an  efficiency  of  12.4%,  a  40  m2  surface  of  PV 
array  is  considered  to  deliver  the  needed  power. 

From  midnight  to  dawn  (0:00-6:30),  the  fuel  cell  is  used  to  cover 
the  power  demand.  When  the  power  delivered  by  the  PV  array 
reaches  the  power  demand  (6:30),  the  PEMFC  generator  is  shut 
down  and  the  PV  generator  takes  charge  of  the  power  requirement. 
During  the  daylight  period,  the  excess  of  power  produced  by  the  PV 
array  is  used  to  produce  hydrogen  from  water  electrolysis.  Hydro¬ 
gen  is  stored  in  a  metal-hydride  tank  to  be  used  for  feeding  the 
PEMFC  generator  when  the  PV  generator  cannot  reach  the  power 
demand  (especially,  during  the  night).  In  the  afternoon,  when  the 
power  delivered  by  the  PV  becomes  lower  than  the  power  demand 
(17:30)  the  PEMFC  generator  starts  up  again  to  cover  the  user  load 
(17:30-24:00). 


2.2.1.  Output  voltage  fuel  cell  modeling 

The  behavior  of  a  fuel  cell  is  strongly  nonlinear  and  depends  on 
several  factors  such  as  the  current  density,  the  fuel  cell  temperature 
and  the  partial  pressure  of  reactants.  The  output  voltage  of  a  single 
fuel  cell  {Vcen)  depends  on  the  thermodynamic  potential  ( ENernst), 
the  ohmic  drop  (p0hm)  and  the  activation  overvoltage  {pact)-  In  the 
present  approach,  the  concentration  overvoltage  links  to  gas  access 
control  is  already  incorporated  in  activation  overvoltage. 

Vcell  —  ENernst  +  Pact  +  Pohm  (1 ) 

The  thermodynamic  potential,  also  called  Nernst  voltage  is 
obtained  at  the  open-circuit  condition  under  thermodynamic  bal¬ 
ance  [15] 

RT 

ENernst  =  1.229-8.5  x  W~4(T  -  298.15)  +  —  In  [Ph2(Po2)°'5]  (2) 

where  T  is  the  fuel  cell  temperature,  Ph2  the  partial  pressure  of 
hydrogen  and  Pq2  the  partial  pressure  of  oxygen.  The  ohmic  over¬ 
voltage  is  written  as  a  function  of  the  temperature  of  the  fuel  cell, 
the  cell  electrode  area  (AfC),  the  membrane  thickness  (Zm)  and  the 
current  (/)  [14]. 

181 .6/[l  +  0.03(I/Afc)  +  0. 062(7/303  )2(//7lfc)2'5](m 
Pohm  —  A/c[11.866-3(//A/c)]exp[4.18((7-303)/7)]  [  ) 


—  Photovoltaic  array  power - Consumption  demand - Fuel  cell  power 


Fig.  2.  Electric  power  produced  by  the  PV  arrays  in  summer  (— )  and  a  hypothetical  charge  profile  for  a  residential  household  in  La  Reunion  Island. 
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The  effects  of  double  layer  capacitance  charging  at  the 
electrode-electrolyte  interface  are  described  by  considering 
the  activation  overvoltage  dynamic  [13]. 

dVact  _  I  Vact  ,as 

~dt~  R^Qi  1  j 

where  Cdi  is  the  double  layer  capacitance,  Ract  =  -  rjactU  the  acti¬ 
vation  resistance.  rjact  the  activation  overvoltage  is  empirically 
modeled  using  the  following  expression  [14]: 

rjact  =  ?i+  &T  +  ln(Co2 )  +  £4rin(J)  (5) 

where 

£2  =  0.00286  +  0.0002  ln(A/c)  +  4.3  x  10“5ln(CH2)  (6) 

§1,  £3  and  §4  are  constant  parametric  coefficients  given  in 
Table  1.  The  oxygen  and  hydrogen  concentrations  at  the  elec¬ 
trode/membrane  interfaces  are  given  using  a  Henry’s  law  equation 
[15]. 

(498  \ 

~^r)  (7) 

Ch2  =  9.174  x  10-7Ph2  exp  (^7)  (8) 

Finally,  the  output  voltage  of  a  stack  (Vstack)  constitutes  of  NceU  cells 
in  series  are  given  by 

^ stack  —  ^cell^cell  (9) 


The  total  power  delivered  is  directly  related  to  the  reacted  hydrogen 
molar  flow  rate  and  the  enthalpy  of  reaction. 

Etot  —  rhy{2reactedA.H  (17) 

The  power  consumed  by  the  electrical  load  and  the  heat  flux  dissi¬ 
pated  at  the  stack  surface  are  given  by 

Eelec  —  V stack.1  (18) 

Efoss  =  hambS(T  —  Tamjj )  (19) 


where  Tamb  is  the  ambient  temperature,  hamb  the  convective  heat 
transfer  coefficient  and  5  the  external  stack  surface.  Finally,  the  heat 
flux  removed  by  the  coolant  is  obtained  using  a  nonlinear  model  of 
the  heat  exchanger  system  [16]. 


ccoo/ 


—  (^cond  +  hconvl) 


(T  —  Tcin)  —  (T  —  Tc,out) 
ln((T-Tc>ftl)/(T-TCioUt)) 


(20) 


where  TCtin  and  Tc,0ut  are  the  input  and  output  water  temperature, 
respectively.  hcond  and  hconv  are  parameters  which  characterize  the 
conduction  and  convection  properties  of  the  heat  exchangers.  Tc,0ut 
is  given  as  a  function  of  the  water  mass  flow  rate  ( ificw ),  the  water 
heat  capacity  (Cpw),  the  water  volume  (Vw),  the  water  density  (pw) 
and  the  overall  heat  transfer  coefficient  (Ua). 


^  \/  r  die, out  _  ^  r  (t  t  \  >  n  It  Tcjn+Tc,out  \ 
Pw V w Cpw — ^ —  —  UicvpLpvv^  lcjn  —  lc,out)  +  c^a  I  1  ^ -  I 


(21) 


2.2.2.  Mass  balance 

According  to  the  ideal  gas  law  and  the  conservation  rule,  the 
variation  of  the  partial  pressure  of  each  gas  depends  on  the  gas 
inlet  flow  rate,  the  gas  consumption  and  the  gas  outlet  flow  rate. 
Thus,  the  dynamic  of  the  partial  pressure  of  hydrogen  (Ph2  )  at  the 
anode  and  of  the  partial  pressure  of  oxygen  (Pq2  )  at  the  cathode  are 
given  by 

^  —  rhy{2  in  —  lflH2  QUt  —  rnn2 , reacted  (19) 

^  —  rh.Q2  in  —  rno2,out  ~  dlQ2  reacted  (H) 

where  mH2  in  and  rh02  in  are  the  inlet  molar  flow  rate  of  hydro¬ 
gen  and  oxygen,  respectively,  mU2  reacted  and  rno2, reacted  the  reacted 
molar  flow  rate  of  hydrogen  and  oxygen,  respectively.  Hydrogen 
and  oxygen  molar  outlet  flow  rates  can  be  estimated  by  the  dif¬ 
ference  between  partial  pressure  inside  the  stack  and  downstream 
pressure  using  the  following  equations  [13]: 

rriu2,out  —  ^an(F>H2  —  Eh2)  (12) 

mo2,out  =  kca(Po2-P°ou2t)  (13) 

The  excess  hydrogen  is  recirculated  to  the  anode.  The  reacted  gas 
flow  rate  is  calculated  by  the  Faraday’s  law  (Eqs.  (14)  and  (15)). 

^H2 , reacted  —  ^cell  2p  ( 1 

^02 ,  reacted  =  N cell  (15) 

2.2.3.  Energy  balance 

The  energy  balance  of  a  PEMFC  is  a  function  of  the  total  power 
{Etot)  from  the  electrochemical  reaction,  the  electrical  power  con¬ 
sumed  by  the  electrical  load  (£e/ec),  the  heat  flux  both  dissipated  to 
the  surroundings  (Pf0ss)  and  removed  by  the  heat  exchanger  (Pcoo/) 

—  Etot  —  Eeiec  —  Eioss  —  Ec00i  (16) 


d  (  ^H2  Van 
dt  ^  RT 

d  (  P02  Vca 
dt  y~RT~ 


All  parameters  values  and  description  are  gathered  in  Table  1. 

3.  Control  design 

In  this  study,  the  PEMFC  generator  is  a  part  of  the  stand-alone 
hybrid  system  and  used  to  cover  the  energy  demand  when  the  PV 
generator  is  deficient.  Consequently,  the  first  controlled  variable  is 
the  output  power  of  the  PEMFC  system.  To  ensure  the  reliability 
of  the  stand-alone  hybrid  system,  the  output  power  of  the  PEMFC 
generator  has  to  track  changes  in  setpoint  with  accuracy  and  short 
time  response  even  when  extreme  load  changes  are  imposed.  To 
reach  these  control  objectives,  the  choice  of  the  manipulated  vari¬ 
able  and  the  controller  are  crucial.  In  this  study,  the  control  of  the 
output  power  is  based  on  a  nonlinear  model-based  predictive  con¬ 
trol  (NMPC)  with  the  molar  flow  rate  of  hydrogen  as  manipulated 
variable.  Several  authors  report  the  robustness  and  efficiency  of  the 
NMPC  when  dealing  with  fuel  cell  control  [9,10,17]. 

It  is  well-known  that  the  performance  and  the  reliability  of 
the  PEMFC  system  strongly  depend  on  its  operating  temperature. 
Therefore,  to  ensure  optimal  operating  condition  and  avoid  any 
risk  of  damage  due  to  high  temperature,  the  control  of  the  PEMFC 
temperature  is  mandatory.  In  this  study,  the  control  of  this  sec¬ 
ond  controlled  variable  is  based  on  a  global  linearizing  control 
(GLC)  approach  using  the  mass  flow  rate  of  water  as  manipulated 
variable.  Differential  geometry  has  been  accepted  as  a  relevant 
tool  for  the  design  of  nonlinear  controllers  and  has  proved  its 
efficiency  for  numerous  nonlinear  chemical  processes.  Most  of 
approaches  concerned  with  this  tool  consist  in  a  feedback  lineariza¬ 
tion,  either  input-output  or  input-state  [18-21].  Regarding  the  GLC, 
an  input-output  linearizing  technique,  the  nonlinear  control  law  is 
based  directly  on  the  nonlinear  model  of  the  PEMFC  system.  Accord¬ 
ing  to  on-line  objective  applications,  GLC  appears  as  a  suitable 
alternative  to  the  widely  used  NMPC  to  reduce  the  computation 
time  of  the  control  variable. 

In  short,  the  control  of  the  output  power  is  performed  using  an 
NMPC  strategy  whereas  the  control  of  the  PEMFC  temperature  is 
based  on  a  GLC  algorithm.  The  benefit  of  the  proposed  control  strat¬ 
egy  (NMPC  +  GLC)  in  terms  of  computation  time  efficiency  is  studied 


148 


C.  Damour  et  al.  /  Journal  of  Power  Sources  206  (2012)  144-152 


Table  1 

Parameter  values  for  the  Ballard  MK5-E  PEMFC  5  kW  system  [13]. 


Symbol 

Description 

Value 

Unit 

Ncell 

Number  of  cell 

35 

Afc 

Cell  electrode  area 

232 

cm2 

h 

Parametric  coefficient 

-0.948 

$3 

Parametric  coefficient 

7.6  xlO-5 

£4 

Parametric  coefficient 

-1.93  x  10-4 

Q 1 

Double  layer  capacitance 

8.12 

F 

Van 

Anode  volume 

0.005 

m3 

Vca 

Cathode  volume 

0.01 

m3 

kan 

Flow  constant  at  the  anode 

0.065 

mol  s-1  atm-1 

kca 

Flow  constant  at  the  cathode 

0.065 

mol  s-1  atm-1 

pin 

h2 

Hydrogen  pressure  at  the  inlet 

3 

atm 

pout 

°2 

Oxygen  pressure  at  the  outlet 

3 

atm 

™02,in 

Oxygen  inlet  flow  rate 

2 

mol  s-1 

lm 

Membrane  thickness 

178  xl0“4 

cm 

hambS 

Stack  heat  transfer  coefficient 

17 

WK"1 

Tamb 

Ambient  temperature 

298 

K 

Tc,in 

Cooling  water  inlet  temperature 

298 

K 

AH 

Hydrogen  consumption  enthalpy 

285.5 

1<J  mol-1 

ct 

Thermal  capacitance 

17.9 

kJK-1 

hcond 

Heat  exchanger  conductive  index 

35.55 

WK”1 

hconv 

Heat  exchanger  convective  index 

0.025 

WK”1  A”1 

Cpw 

Heat  capacity  of  water 

4.184 

kjkg-1  K_1 

Pw 

Water  density 

1000 

kg  m-3 

vw 

Water  volume 

2.5  x  10-3 

m3 

Ua 

Overall  heat  transfer  coefficient 

241 

WK-1 

in  Section  4.  With  this  aim,  a  comparison  between  the  proposed 
strategy  and  a  classical  multivariable  NMPC  strategy  is  performed. 

3A.  NMPC  of  output  fuel  cell  power 

The  NMPC  is  based  on  a  prediction  of  the  future  output  of 
the  system  to  calculate  the  current  control  action.  In  practice,  the 
current  control  action  is  obtained  by  solving  on-line  an  optimiza¬ 
tion  problem.  The  aim  of  the  optimization  problem  is  to  find  the 
optimal  solution  of  a  nonlinear  cost  function  that  minimizing  the 
mean  squared  difference  between  predicted  outputs  and  target  val¬ 
ues.  This  strategy  offers  a  great  flexibility  and  robustness.  In  fact, 
the  optimization  problem  can  handle  penalty  for  variations  in  the 
control  variable  (e.g.  to  avoid  excessive  variation),  degree  of  impor¬ 
tance  of  the  setpoint  tracking  accuracy,  the  error  modeling,  etc.  In 
this  study,  the  optimization  problem  can  be  rewritten  as 

min  J(k)  =  min 

u(k/k),u(k+l/k),...,u(k+Nu/k)  u(k/k),u(k+l/k),...,u(k+Nu/k) 

(Ny  Nu  \ 

^y(e(k  +  i/k))2  +  ^pAu2(k  +  i/k )  (22) 

>'=n  i  «'= i  ) 

subject  to  constraints  on  the  manipulated  variable 

Umm<u(k  +  i/k)<umax  for  i  =  l,2,  ...,NU  (23) 

where  u(/<),  the  control  variable,  is  the  molar  flow  rate  of  hydrogen 
at  time  /<,  Ny  and  Nu  the  predictive  and  control  horizon,  respec¬ 
tively,  y  and  /3  two  weighting  parameters.  The  prediction  horizon 
corresponds  to  the  future  time  interval  used  to  compute  the  pro¬ 
cess  output  predictions  with  the  PEMFC  model.  The  control  horizon 
corresponds  to  the  time  interval  when  present  and  future  control 
actions  are  computed.  The  optimization  parameter  Ni  determines, 
together  with  the  prediction  horizon  the  coincidence  horizon.  The 
trajectory  tracking  error  and  the  changes  in  the  control  variable  are 
formulated  by 

e(/<  +  i/k)  =  ysr(k  +  i/k )  -  (ym(fc  +  i/k)  +  y(/<)  -  ym{k))  (24) 

A u(k  +  i/k)  =  u(k  +  i/k)  -  u(k  +  i  -  1  //<)  (25) 


where  ysr  is  the  target  value  and  ym  the  model  output.  ysr  allows  to 
take  into  account  the  dynamic  of  the  system,  it  gives  the  trajectory 
to  follow  to  make  the  process  output  reach  the  setpoint  w.  Fig.  3 
illustrates  the  overall  NMPC  principle. 


3.2.  GLC  of  fuel  cell  temperature 

The  objective  of  input-output  linearization  is  to  obtain  a  non¬ 
linear  control  law  directly  based  on  the  nonlinear  model  of  the 
system.  To  understand  the  GLC  approach  proposed  in  this  paper,  a 
brief  recall  about  the  concept  of  differential  geometry  is  essential. 
For  brevity,  it  will  be  limited  to  single-input  single-output  (SISO) 
system.  However,  the  idea  can  be  straightforwardly  extended  to  a 
MIMO  system.  Considering  the  nonlinear  SISO  system  given  by 

*(t)  =  /(x)  +  g(x)u  (26) 

y(x)  =  WO)  (27) 


Fig.  3.  NMPC  principle. 
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Fig.  4.  Overall  control  scheme  of  the  PEMFC  system,  where  yiset  and  y2set  are  the  output  power  setpoint  and  the  temperature  stack  setpoint,  respectively. 


where  x  e  Rn  is  the  state  vector,  u  the  input  vector  and  y  the  output 
vector.  An  input-output  relation  can  be  deduced 

yr(t)  =  Lfh(x(t))  +  Lg(Lg~J  h(x(t)))u(t)  (28) 

where  yr(t)  =  dry/dtr,  LJ/i  and  Lrgh  are  the  rth  Lie  derivatives  of  h(x(t)) 
along/and  g.  r  <  n  is  the  relative  degree  of  the  system,  i.e.  the  num¬ 
ber  of  times  the  output  has  to  be  differentiated  with  respect  to  time 
before  one  input  u  appears  explicitly  in  the  resulting  equations. 

Once  this  relation  is  established,  a  closed-loop  control  can 
be  designed.  In  most  process  control  applications,  to  reduce  the 
effect  of  unmeasured  disturbances  and  plant/model  mismatch,  the 
control  law  includes  an  integral  term  [18,22,23].  Consequently, 
considering  a  setpoint  to  track,  denoted  ysp,  the  tracking  problem 
is  given  by 


yr(t)  =  ysrP  + 


r-1 

^k;9f«(t)  +  ao 

i=0 


/' 


m(t)dz 


(29) 


where  9f(  t)  =ysp(  t)  -  y(  t)  is  the  pursuit  error,  St®  =  d1 9t/dt* ,  k0,  fci . 

kr_ i  and  «o  the  controller  tuning  parameters.  Finally,  the  control 
law  is  deduced  from  Eqs.  (28)  and  (29) 


u(t)  = 


yrsP  +  + a0  j*  mdr  -  lj/imo) 

LgiVf'wm 


(30) 


In  this  study,  according  to  the  dynamic  nonlinear  model  of  the 
PEMFC  system,  the  GLC  approach  is  used  to  control  the  temperature 
of  the  PEMFC  with  the  water  mass  flow  rate  as  control  variable. 


3.3.  Control  of  the  PEMFC  system 

The  control  of  the  PEMFC  system  is  based  on  a  NMPC  of  the  out¬ 
put  power  (yi )  using  the  hydrogen  molar  flow  rate  ( U\ )  as  control 
variable  and  on  a  GLC  strategy  to  regulate  the  PEMFC  tempera¬ 
ture  (y2)  by  manipulated  the  water  mass  flow  rate  (u2).  Widely 
known  and  used  for  its  convergence  and  robustness  properties,  the 
Levenberg-Marquard’s  algorithm  is  used  here  to  solve  on-line  the 
optimization  problem  of  the  NMPC  strategy.  The  overall  control 
scheme  of  the  MIMO  PEMFC  system  is  illustrated  in  Fig.  4. 


Due  to  physical  constraints,  the  water  mass  flow  rate  range  from 
0.3  to  0.12  kg  s-1  and  the  hydrogen  molar  flow  rate  range  from  0  to 
1  mols-1. 

4.  Simulation  results 

In  order  to  evaluate  the  performance  of  the  proposed  control 
scheme,  a  5  kW  PEMFC  plant  system  is  simulated  in  the  Matlab™ 
environment.  It  should  be  pointed  out  that  an  operating  fuel  cell 
plant  requires  the  regulation  of  other  variables,  to  avoid  for  instance 
oxygen  starvation  or  membrane  flooding  or  dehydration.  In  the  fol¬ 
lowing  simulations,  it  is  assumed  these  variables  have  already  been 
well  controlled. 

4.1.  Control  performance 

The  performance  of  the  proposed  control  strategy,  in  terms  of 
setpoint  tracking,  disturbances  rejection  and  robustness  against 
model  mismatch  in  presence  of  noise  measurement,  is  studied  via 
three  control  scenarios.  The  first  scenario  evaluates  the  perfor¬ 
mance  of  the  proposed  strategy  when  large  step  changes  of  load 
power  demand  appear.  The  second  one  allows  us  to  study  the 
proposed  strategy  when  a  step  temperature  is  set.  Finally,  to  eval¬ 
uate  the  robustness  of  the  controller,  the  third  scenario  considers  a 
model  mismatch  in  parameters.  In  this  purpose,  the  flow  constant 
at  the  anode  ( kan )  and  the  overall  heat  transfer  coefficient  (U)  are 
randomly  increased  and  decreased  from  their  nominal  values  by 
20%  every  10  min.  Moreover,  to  be  closer  to  the  experimental  con¬ 
ditions  a  white  Gaussian  noise  with  a  signal  to  noise  ratio  of  1%  is 
added  to  the  process  outputs  simulator. 

Fig.  5  shows  results  for  the  first  scenario.  A  good  performance  in 
PEMFC  output  power  setpoint  tracking  and  a  satisfactory  regulation 
of  the  stack  temperature  in  presence  of  power  disturbances  are 
established.  Moreover,  the  dynamic  of  both  manipulated  variables 
is  entirely  suitable.  This  first  scenario  highlights  the  requirement 
of  an  appropriate  control  of  the  stack  temperature  to  avoid  any 
damage  of  the  PEMFC. 

Fig.  6  shows  a  correct  running  of  the  proposed  strategy  in 
temperature  setpoint  tracking,  considering  a  step  change  of  fif¬ 
teen  degrees,  and  an  excellent  regulation  of  the  output  power  in 
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Fig.  5.  First  control  scenario:  (a)  PEMFC  output  power;  (b)  stack  temperature;  (c)  hydrogen  molar  flow  rate;  (d)  water  mass  flow  rate. 
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Fig.  6.  Second  control  scenario:  (a)  PEMFC  output  power;  (b)  stack  temperature;  (c)  hydrogen  molar  flow  rate;  (d)  water  mass  flow  rate. 
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Fig.  7.  Third  control  scenario:  (a)  PEMFC  output  power;  (b)  stack  temperature;  (c)  hydrogen  molar  flow  rate;  (d)  water  mass  flow  rate. 


presence  of  temperature  disturbances.  The  saturation  observed  in 
Fig.  6d  between  28  and  33  min  is  due  to  previously  mentioned 
physical  constraints  on  water  mass  flow  rate.  This  second  scenario 
confirms  that  the  performance  of  the  PEMFC  strongly  depends  on  its 
operating  temperature.  In  one  hand,  too  lower  temperature  slows 
down  mass  and  charge  transfer  reducing  the  PEMFC  performance. 
In  the  other  hand,  too  higher  temperature  could  damage  the  PEMFC 
and  impacts  its  performance. 

Fig.  7  demonstrates  the  control  effectiveness  of  the  multivari¬ 
able  proposed  control  strategy  under  several  step  changes  of  high 
load  power  demand  in  presence  of  high  level  of  noise  measure¬ 
ment  and  plant/model  mismatch.  The  excellent  accuracy  in  power 
setpoint  tracking  and  the  satisfactory  regulation  of  the  stack  tem¬ 
perature  emphasize  the  robustness  of  the  proposed  controller 
against  unavoidable  plant/model  mismatch  and  noise  measure¬ 
ment. 

Flowever,  whereas  the  manipulated  variable  computed  by  the 
NMPC  controller  is  smooth,  it  can  be  pointed  out  that  the  water 
flow  rate  calculated  using  GLC  is  affected  by  the  noise.  In  the  case 
of  GLC,  since  the  control  law  is  directly  deduced  from  the  model 
process,  the  manipulated  variable  is  affected  by  the  noise  measure¬ 
ment.  In  practice,  this  limitation  is  easily  overcome  either  filtering 


on-line  the  measurements  or  directly  the  control  action.  A  widely 
used  heuristic  rule  consists  of  applying  a  first-order  filter  to  the 
control  action  to  reduce  high  frequency  oscillation.  In  this  manner, 
the  filtered  control  action  i/dtered  is  obtained  using  the  relation  (31): 

ufered(t)  =  ca4ltered(t  -  1)  +  (1  -  a)ufc{t)  (31) 

where  u^LC(t)  is  the  control  action  calculated  by  the  GLC  approach 
and  a  the  filter  coefficient  between  0  and  1. 

Flowever,  when  dealing  with  optimal  control,  this  approach 
leads  to  a  sub-optimal  solution  of  the  control  action,  which  in  some 
cases  could  involve  a  significant  decrease  of  the  setpoint  tracking 
accuracy. 

In  the  present  case,  an  appropriate  setting  of  the  filter  param¬ 
eter  allows  to  smooth  the  control  action  without  decreasing  the 
setpoint  tracking  performance.  Considering  a  filtering  parameter  a 
equal  to  0.5,  Fig.  8  illustrates  the  benefits  of  this  heuristic  rule  to  sig¬ 
nificantly  reduce  the  oscillation  of  the  control  variable  in  presence 
of  noise  measurement  and  plant/model  mismatch  (third  scenario 
configuration). 


a 


Time  /  h 


Fig.  8.  Filtering  approach  performance:  (a)  stack  temperature;  (b)  water  mass  flow  rate. 
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Table  2 

Comparison  between  proposed  control  strategy  and  classical  NMPC  strategy  in 
terms  of  accuracy  and  computational  time. 


Classical  strategy 

Proposed  strategy 

Profile  1 

RMSE  for  power 

0.468 

0.329 

RMSE  for  temperature 

0.253 

0.176 

MCT  (s) 

20.028 

4.76 

Profile  2 

RMSE  for  power 

0.812 

0.687 

RMSE  for  temperature 

0.274 

0.218 

MCT  (s) 

18.479 

5.080 

Profile  3 

RMSE  for  power 

0.787 

0.712 

RMSE  for  temperature 

0.365 

0.318 

MCT  (s) 

21.374 

5.728 
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5.  Conclusions 

In  this  work,  an  innovative  model-based  control  approach  of 
a  proton  exchange  membrane  fuel  cell  is  proposed.  In  this  aim,  a 
MIMO  dynamic  nonlinear  model  of  a  PEMFC,  dedicated  to  nonlinear 
model-based  control  approaches,  is  developed  and  implemented  in 
Matlab™.  To  control  the  output  power  and  the  stack  temperature, 
a  nonlinear  model-based  predictive  control  strategy  (NMPC)  and 
a  global  linearizing  control  (GLC)  approach  are  designed,  directly 
based  on  this  model. 

The  simulation  results,  for  a  5  kW  PEMFC  power  plant,  demon¬ 
strated  that  the  proposed  control  strategy  led  to  good  performance 
in  setpoint  tracking,  disturbances  rejection  and  robustness  against 
plant/model  mismatch  and  noise.  Regarding  on-line  application 
objective,  the  proposed  strategy  appears  to  be  four  times  faster  than 
classical  multivariable  NMPC  strategies.  It  is  well-known  that  the 
cell  performance  changes  along  its  life,  the  model  parameters  have 
to  be  adjusted  to  reduce  plant/model  mismatch.  Therefore,  semi- 
empirical  PEMFC  model  used  by  this  control  approach  makes  easy 
the  tuning  of  the  model  parameters  in  order  to  reduce  plant/model 
mismatch  during  lifetime  of  cell.  The  tuning  of  the  model  param¬ 
eter  could  be  performed  off-line  through  experiments,  or  on-line 
using  an  optimization  algorithm. 

This  proposed  model-based  control  approach  appears  suitable 
to  be  included  in  the  overall  control  scheme  of  the  stand-alone 
hybrid  power  generation  system  used  to  cover  the  daily  power 
consumption  of  a  typical  household  in  subtropical  island  (lie  de 
La  Reunion:  21°  6'  S/55°  36r  E). 
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In  order  to  highlight  the  benefit  due  to  reduction  of  com¬ 
putational  time,  classical  multivariable  NMPC  strategy  has  been 
compared  to  the  proposed  control  strategy  (NMPC  +  GLC).  There¬ 
fore,  a  classical  multivariable  NMPC  strategy  is  implemented  in 
Matlab®  environment  considering  the  water  mass  flow  rate  and 
the  hydrogen  molar  flow  rate  as  manipulated  variables,  and  the  out¬ 
put  power  and  the  stack  temperature  as  controlled  variables.  The 
performance  of  each  strategy,  in  terms  of  computational  time  effi¬ 
ciency,  is  evaluated  on  three  typical  power  demand  profiles  (Fig.  9), 
considering  a  constant  setpoint  of  70  °C  for  the  stack  temperature. 

Table  2  summarizes  the  computational  time  efficiency  and  the 
setpoint  tracking  accuracy  of  each  strategy.  This  study  proposes 
simulations,  using  Matlab®  software  and  running  with  2-GHz  and 
32-bit  computer.  The  root  mean  square  error  (RMSE)  criterion  is 
used  to  evaluate  the  performance  in  terms  of  setpoint  tracking 
accuracy,  whereas  the  computational  time  efficiency  is  investigated 
using  the  mean  computational  time  (MCT): 


RMSE  = 


MCT  = 


E"=lCTi 


(32) 

(33) 


where  n  is  the  number  of  data  and  CT the  time  required  to  solve 
the  control  problem  at  time  t=i. 

It  can  be  noticed  that  for  similar  accuracies  of  RMSE,  the  com¬ 
bination  of  NMPC  and  GLC  algorithms  is  four  times  faster  than 
a  classical  multivariable  NMPC  strategy.  Regarding  on-line  appli¬ 
cation  objectives,  the  proposed  control  strategy  appears  to  be  a 
suitable  alternative  to  the  classical  multivariable  NMPC  controller. 
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